Meissner effect in diffusive normal metal / d-wscve superconductor 

junctions 

Takehito Yokoyama'^'^, Yukio Tanaka^, Alexander Golubov^ Jun-ichiro Inoue^, 

Yasuhiro Asano ^ , 

'^Department of Applied Physics, Nagoya University, Nagoya 464-8603, Japan. 
CREST Japan Science and Technology Corporation (JST), Nagoya 464-8603, Japan. 
^ Faculty of Science and Technology, University of Twente, 7500 AE, Enschede, The Netherlands. 
Department of Applied Physics, Hokkaido University, Sapporo 060-8628, Japan. 



Abstract 

The Meissner effect in diffusive normal metal / insulator / d-wave superconductor junctions is studied theoretically 
in the framework of the Usadel equation under the generalized boundary condition. The effect of midgap Andreev 
resonant states (MARS) formed at the interface of d-wave superconductor is taken into account. It is shown that 
the formation of MARS suppresses the susceptibility of the diffusive normal metal. 
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1. Introduction 

In diffusive normal metal / superconductor 
(DN/S) junctions, the DN acquires induced su- 
perconductivity, i.e. Cooper pairs penetrate into 
the DN. This proximity effect has been studied 
since the BCS theory was established. The prox- 
imity induced Meissner demagnetization in DN/S 
junctions was measured experimentally by Oda 
and Nagano[l] and Mota et al.[2]. It has T'^/^ 
dependence in the dirty limit. The quasiclassical 
Green's function theory was used earlier to study 
the Meissner effect in proximity structures. 
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The quasiclassical Green's function theory was 
developed by Eilenberger [3] and was generalized 
by Eliashberg [4], Larkin and Ovchinnikov [5] in 
order to study the noncquilibrium state. This the- 
ory was applied by Zaikin[6] and Kieselmann[7] 
to studing the Meissner effect in DN/S junctions. 
Narikiyo and Fukuyama [8] calculated the Meiss- 
ner screening length in a semi-infinite system con- 
taining an Anderson impurity. Higashitani and Na- 
gai studied the Meissner effect in the clean limit [9] . 
Belzig et al. [10,11] have considered more realistic 
systems by assuming a perfectly transparent N/S 
interface. Up to now the boundary conditions de- 
rived by Kupriyanov and Lukichev (KL) [12] were 
widely used to study proximity effect in DN/S 
structures. 

A more general boundary conditions was de- 
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rived by Nazarov [13] based on the Keldysli-Nambu 
Green's function formalism [14] within the frame- 
work of the Landauer-Biittiker scattering formal- 
ism. The merit of this boundary condition is that 
the BTK theory [15] is reproduced in the ballistic 
limit while in the diffusive limit with a low trans- 
missivity of the interface, the KL boundary con- 
dition is reproduced. Although almost all previ- 
ous papers on Meissner effect in mesoscopic NS 
junctions are either based on the KL boundary 
conditions or on the BTK model, in the actual 
junctions, the transparency of the junction is not 
always small and the impurity scattering in the 
DN cannot be neglected. Tanaka et al.[16] and 
Yokoyama et al.[17] calculated tunneling conduc- 
tance by using the Nazarov's boundary condition. 

It is well known in d-wave superconductors that 
the midgap Andreev resonant states (MARS) are 
formed at the interface of d-wave superconduc- 
tor. The MARS crucially influence various physi- 
cal quantities [18]. One of the authors (Y.T.) re- 
cently generalized the boundary condition of the 
Keldysh-Nambu Green's function formalism to un- 
conventional superconductor junctions [19,20]. It is 
revealed that in DN/d-wave superconductor junc- 
tions the proximity effect and the MARS strongly 
compete with each other [19], while they coex- 
ist in DN/triplet superconductor junctions. The 
newly obtained boundary conditions expressed in 
the Keldysh-Nambu Green's function are useful for 
the calculation of various physical quantities. The 
timely problem is to study theoretically the Meiss- 
ner effect in DN / d-wave S junctions using the new 
boundary conditions [19]. In the present paper, we 
calculate the susceptibility of the DN layer in DN/ 
d-wave S junctions for various junction parameters 
such as the height of the insulating barrier at the 
interface and the angle between the normal to the 
interface and the crystal axis of a d-wave supercon- 
ductor. 

The organization of the paper is as follows. In 
section 2, we will provide the derivation of the ex- 
pression for the susceptibility of the DN. In section 
3, the results of calculation are presented for vari- 
ous types of junction. In section 4, the summary of 
the obtained results is given. In the present paper 
we set c — kB ~ h ^ 1. 



2. Formulation 

In this section, we introduce the model and the 
formalism. We consider a junction consisting of 
vacuum (VAC) and superconducting reservoirs 
connected by a quasi-one-dimensional diffusive 
conductor (DN) with a length L much larger than 
the mean free path. We assume that the interface 
between the DN conductor and the S electrode at 
X = L has a resistance the DN/VAC interface 
at a:; = is specular, and we apply the general- 
ized boundary conditions by Tanaka [19] to treat 
the interface between DN and S. A weak external 
magnetic field H is applied in z-direction (see Fig. 
1). The vector potential can be chosen to have 
only the y component which depends on x. 
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Fig. 1. Schematic illustration of the model. 

We describe the insulating barrier between DN 
and S by using the (5-function (i.e., U (x) — H6{x — 
L)), which provides the transparency of the junc- 
tion Trn = 4cos^ 0/(4cos^ (/) -I- Z^), whcrc Z = 
2H/vf is a dimensionless constant, (j) is the injec- 
tion angle of a quasiparticle measured from the in- 
terface normal to the junction and vp is Fermi ve- 
locity. 

In the following, we solve the Usadel equations 
[21] with using the standard 6'-parameterization. 
The parameter 9{x) is a measure of the proximity 
effect in DN and obey the following equation 
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D—9{x) - 2ujnsm[9{x)] = 0, 
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where D and a;„ denote the diffusion constant and 
the Matsubara frequency, respectively. The bound- 
ary condition for 0{x) at the DN/S interface is 
given in Ref. [19]. The interface resistance Rb is 
given by 



Rb — Rq- 



PJ' d<i>T{<j>) cos 



with T{(j)) = 4cos2(/)/(4cos2 
-1 



(2) 



^2). Here i?o is 



Sharvin resistance = e^fc|,S'c/(47r^), where fc^ 
is the Fermi wave number and Sc is the constriction 
area. The current distribution is given by 
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Hx)A{x), (3) 



where A{x), N{0) and T denote the vector poten- 
tial, the density of states at the Fermi energy and 
the temperature of the system respectively. The 
Maxwell equation reads 



A (x) = -Anj (x) . 



(4) 



The boundary conditions for A(x) are given by 
d 



dx 



A{0)^H, A{L)^0, 



(5) 



where we have neglected the penetration of mag- 
netic fields into the superconductor by assuming a 
small penetration depth in S. 

Finally we obtain the expression of the suscep- 
tibility. 



Attx = 1 
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The d-wave pair potentials in directional space are 
given by A± = A(T) cos2((?!) =F ct), where A(r) is 
the magnitude of pair potential at a given temper- 
ature T and a denotes the angle between the nor- 
mal to the interface and the crystal axis of a d-wave 
superconductor. 



3. Results 

In the following, we focus on the magnitude of 
the diamagnetic susceptibility x induced by the 



proximity effect. Figs. 2 and 3 show the suscepti- 
bility for Z = 10 and Z = respectively where 
K = IGne^ N (0) . For a = 0, the temperature 
dependencies of — 47rx are not much different. For 
a = 0.1257r, the magnitude of x for Z = 10 is much 
stronger suppressed than that for Z = 0. At the 
same time, we find that the magnitude of x de- 
creases with increasing a. We note in the case of 
a = 0.257r that the susceptibility completely van- 
ishes, (i.e., —4:TTx = 0). This is because the prox- 
imity effect is absent in diffusive metals due to an- 
gular averaging [19]. The absence of the proximity 
effect is a significant feature specific for junctions 
containing unconventional superconductors. 
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Fig. 2. Susceptibility for low transparent junctions with 
Z = 10. 
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Fig. 3. Susceptibility for high transparent junctions with 
Z = 0. 

We also plot the a dependencies of the suscep- 
tibility at T/Tc = 0.01 and T/Tc = 0.1 in Fig. 
4. For all cases, x is a decreasing function of a. 
At T/Tc = 0.01, the magnitude of x for Z ^ 10 
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rapidly decreases with the increase of a. The re- 
sults imply that the MARS suppresses the prox- 
imity effect in low transparent junctions and low 
temperatures. 
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Fig. 4. a dependences of the susceptibility at 
T/Tc = 0.01{upper panel) and T/Tq = 0.1{lower panel). 



4. Conclusions 

In the present paper, we have calculated the 
induced Meissner effect by the proximity effect in 
DN region of DN/d-wave superconductor junc- 
tions. We have solved the Usadel equation under 
a general boundary condition [19] in which the 
formation of the MARS is fully taken into ac- 
count [18]. The magnitude of x decreases with the 
increase of a up to 0.257r. At a = 0.257r, where all 
quasiparticles feel MARS, the x becomes zero. It 
might be interesting to check experimentally such 
an anomalous proximity effect in DN. Another fu- 
ture problem is a similar calculation of the induced 
Meissner effect with a p-wave triplet superconduc- 
tor instead of a d-wave one, since dramatic new 
phenomena are recently predicted in DN/ triplet 
junctions [20]. 
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